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Abstract 

We study the backreaction problem in a mechanism of magnetogenesis from inflation. In usual 
analysis, it has been assumed that the backreaction due to electromagnetic fields spoils infiation 
once it becomes important. However, there exists no justification for this assumption. Hence, we 
analyze magnetogenesis from inflation by taking into account the backreaction. On the contrary 
to the naive expectation, we show that inflation still continues even after the backreaction begins 
to work. Nevertheless, it turns out that creation of primordial magnetic flelds is significantly 
suppressed due to the backreaction. 
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I. INTRODUCTION 



The origin of cosmological magnetic fields still remains an enigma in modern precision 
cosmology [ll, l2l]- In particular, how to make magnetic fields on Mpc scales would be a big 
challenge. As an attractive possibility, creation of primordial magnetic fields from infla- 
tion has been discussed s], 0] . In particular, a nonminimal kinetic term of vector fields in 
supergravity can be used to generate primordial cosmological magnetic fields j^. 

Recently, the backreaction issue in this mechanism is raised |6|. There, it is pointed out 
that magnetic fields from inflation may not be a viable mechanism when the backreaction is 
taken into account. The difficulty to create magnetic fields from inflation originates in the 
existence of electric fields whose energy density rapidly grows and soon catches up with the 
energy density of inflaton. Thus, the backreaction begins to work. The crucial assumption 
here is that the backreaction destroys inflation. Hence, creation of magnetic fields cannot 
be compatible with a sufficiently long inflation. However, no one has proved the above 
assumption. Rather, there is a counter example that the backreaction of vector fields does 
not spoil inflation 7|. Hence, it is worth studying magnetogenesis from inflation by taking 
into account the backreaction. 

In this paper, we reconsider generation of magnetic fields during inflation. We take into 
account the backreaction of electromagnetic fields in an inflationary scenario. It turns out 
that the resultant cosmological evolution with backreaction is drastically different from the 
one without backreaction. Indeed, the dynamics of inflaton is quite different and the universe 

Q 

is anisotropically inflating although anisotropy is very small [7|]. In this new background, 
we consider generation of magnetic fields and conclude that cosmological magnetic fields 
could be hardly produced during inflation due to the backreaction although inflation still 
continues in spite of backreaction. 

The organization of the paper is as follows. In section II, we present our model and 
clarify the condition when the backreaction becomes important. In section HI, we review a 
naive perturbative analysis and reveal the necessity of the backreaction. In section IV, we 
take into account the backreaction and show that the backreaction changes the dynamics of 
inflaton significantly. In this new inflationary background, we study generation of magnetic 
fields and calculate the power spectrum of magnetic fields. We find that it is difficult to 
produce magnetic fields from inflation due to backreaction although inflation still continues. 
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The final section is devoted to conclusion. 



II. MODELS 

In this section, we present our models and basic equations for discussing the backreac- 
tion 0]. Here, we envisage the chaotic inflation although it is easy to extend analysis to 
other scenarios. We consider electromagnetic fields in a homogeneous anisotropic universe 
and explain why it is believed that inflation is spoiled when the backreaction of electro- 
magnetic fields becomes important. Then, treating electromagnetic fields perturbatively, we 
clarify in which case the backreaction becomes important. 

Let us consider the following action for the gravitational field, the inflaton field (p and 
the electro-magnetic vector field coupled with 0: 



(1) 



where Mp is the reduced Plack mass, g is the determinant of the metric, R is the Ricci scalar, 
V{(f)) is the inflaton potential, /(0) is the coupling function of the inflaton field to the vector 
one, respectively. The field strength of the vector field is defined by F^^, = d^A^, — d^A^. 

Thanks to the gauge invar iance, we can choose the gauge = 0. Without loss of 
generality, we can take x-axis in the direction of the vector. Hence, we take the homogeneous 
fields of the form = ( 0, Ax{t), 0, ) and cf) = (f){t) . Note that we have assumed the 
direction of the vector field does not change in time, for simplicity. This field configuration 
holds plane symmetry in the plane perpendicular to the vector. Then, we take the metric 
to be 

ds^ = -dt' + e^^W [ e-^-^dx^ + e^'^W {dy' + dz') ] , (2) 

where the cosmic time t is used. Here, a = e" is an isotropic scale factor and a represents 
a deviation from the isotropy. It should be noted that we need to consider an anisotropic 
spacetime from consistency when we treat a background vector field. With above ansatz, 
one obtains the equation of motion for the vector field which is easily solved as 

= /-'(0)e-"-^>^, (3) 

where an overdot denotes a derivative with respect to the cosmic time t and pa denotes a 
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constant of integration. Substituting ([3]) into other equations, we obtain basic equations 

(4) 



• 2 -21 ^ 

a = a + 



3M2 



2" 2 



^--3«^ + ^n0) + ^/-(0)e— (5) 

= -3dd+^/-^(0)e-^'^-^-, (6) 
= -3d0 - + p^/-^(0)/^(0)e-^°-^<^ , (7) 

where denotes the derivative of V with respect to cj). 

From Eq.(jll), we see the effective potential Ves = V + p^/~^e"^°~^'^/2 determines the 
inflaton dynamics. As the second term is coming from the vector contribution, we refer it to 
the energy density of the vector. Let's check if inflation occurs in this model. Using Eqs.(jl]) 
and ([5]), equation for acceleration of the universe is given by 



a + a= -2&^ - 



a my 3M^ 



-.2 



Pa r-2„-4Q-4o- 
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(8) 



It is easy to see that if the energy density of the vector is dominant inflation soon ends. This 
is the reason why it is usually supposed that the backreaction spoils inflation. However, 
as we will see in section IV, the relevance of backreaction does not necessarily imply the 
energetic dominance of the electromagnetic fields. 

To see when the backreaction of electromagnetic fields becomes relevant, we consider 
electromagnetic fields in the isotropic universe. From Eq.(j4]), it is apparent that the fate 
of the electric field depends on the behavior of coupling function By considering the 

critical case /(0) oc e~^" for which the energy density of the vector field remains almost 
constant during inflation, we can determine the functional form of / under the assumption 
of slow-roll inflation. We can use conventional slow-roll equations 

«' = ^n0), Sa<P = -VM. (9) 

Then, we have an equation da/dcf) = a/cj) = —V{(j))/MpV^{(j)) . This can be easily integrated 
as a = — J V/MpV^d(f) . Here, we have absorbed a constant of integration into the definition 
of a. Thus, we obtain 



2 r V 



f = e-'"^ = e"p ""^"^ . (10) 



4 



For the polynomial potential V oc (j)'^, we have / = e'^'^^'^^p . Let us consider more general 
cases by introducing a parameter c in the following form j^: 

/ = e^"^' = oc a-^' , (11) 

where we have used the relation 

gi^'^' = . (12) 

We notice that the energy density of the vector field during inflation would be negligible 
when c < 1 and remains constant when c = 1. While it grows when c > 1 and the growth 
rate of the energy density of vector fields can be calculated as 

^/-2(0)e-^°ocexp(4(c-l)a) . (13) 

Apparently, we need to take into account the backreaction for c > 1. In the next section, 
we consider generation of magnetic fields in the conventional inflationary background driven 
by an inflaton. There, we will recognize that c > 1 is necessary for creating cosmological 
magnetic fields and therefore the backreaction turns out to be important. 



III. MAGNETIC FIELDS FROM INFLATION 



In this section, we review the arguments made in the paper by Demozzi et al. 
Let us consider the Maxwell fields 

1 



S 



d xy/—g 



(14) 



where we assumed a general coupling function (ITTI) . We will consider the isotropic inflation- 
ary background 



ds^ = a^irj) {-dTf + dx^ + dy^ + dz^) 



(15) 



where, for convenience, we switched to the conformal time rj. It is well known that the 
physical degrees are described by the transverse vector which can be written in Fourier 
space as 

d^k 



(16) 



where a denotes degrees of polarization and polarization vectors ef (k) satisfy the relation 
kie1{]s) = and e°'(— k)ef(k) = ^^.p • Assuming the isotropic FRW universe f[T5l) . we can 
deduce the action for physical modes 

^ = ^ E / dvd'kf{<P) [A'^A'^\ - k'AlA'L^] , (17) 

(T=l,2 

where a prime represents a derivative with respect to the conformal time. Let us quantize 
the system by promoting A^ and the conjugate momentum vrj = to operators 

satisfying [A^, tt^,] = i6'^^^6(k — k') , and others are zero. Using creation and annihilation 
operators satisfying a^, a^l = 5'^''5(k — k') , we can expand A^ as 

Al = u^al + n*ai , (18) 

where mode functions have to satisfy equations of motion 

f 

K + 2y< + = (19) 

and the normalization conditions u^u^^ — ulji\^ = i/ P ■ From the correlation function of the 
vector field 

/rlh f rlh 

T;?^''"^''"^''^ = j ^^A{k,v)e'''- , (20) 

we can read off the power spectrum 

Sl{k,v)=^-^. (21) 



Similarly, we can deduce the power spectrum of the magnetic fields Bi = eijkF^ /2 

Sl{k,v) = H^^ (22) 



where we used the formula 

B'Bi = ^FikFik = ^ {diAkdiAk - dkAAAk) . (23) 

The vacuum state |0 > can be specified by the initial condition for the positive frequency 
mode at a sufficiently past time rj 



The above mode function on subhorizon scales connects to the superhorizon solutions 



= C^ + C, I ^, (25) 



where Ci and C2 are constants of integration. For the coupling function, we take 

where Ces is a parameter and a/ is the scale factor at the end of inflation rjf. In the case of 
the conventional slow roll inflation, following (fTT]) . we have Cefr = c. Then, the solution (!25|l 
becomes 

u^ = Ci + C2a^'-^ , (27) 

where we used the relation drj = da/{Hja^) and rescaled constants of integration. Here, the 
Hubble parameter during inflation Hj is assumed to be constant. There are two branches 
where magnetic fields can be created. For a negative c, there exists no backreaction problem 
because the electric fields do not exist in this case. For this case, however, the analysis is 

n 

not reliable due to the strong coupling problem [6|. Hence, we will not consider this case 
hereafter. For a positive c, the second term of (127|) which depends on time is a relevant 
one. By matching and TI7}i at the horizon crossing a^Hj = k, the constant C2 can be 
determined. Thus, we obtain 

where we defined = [a^ / a f)"'^'^ ■ Substituting (l28l) into (!22|) . we obtain magnetic fields at 
the end of inflation ?7j as 

where Ap = af/k is the physical wavelength at the end of inflation corresponding to the 
comoving wavenumber k. To get the flat spectrum, we need c = 3/2. Assuming the GUT 
scale inflation Hi ~ 10~^Mp, we get 5b ~ lO^^^M^ at the end of inflation. Since IG ~ 
10-2°GeV^ this implies 5b ~ lO'^'^G at the end of inflation. If we assume the instantaneous 
reheating and take {MpHiY^'^ as a reheating temperature, we have 



a, T„ ^ ' 
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where we used the temperature Tq ~ 10~^^ GeV at present oq. Taking into account the 
relation 6b oc 1/a^, we obtain 5b(Ap,?7o) ~ {df /aoY6B{^p,r]o) ~ lO^^^G at present. This 
is close to the expected value lO^^G from observations. In that case, however, the energy 
density of the electromagnetic fields Pem exceeds the energy of the inflaton. In fact, from 
the formula for the energy-momentum tensor 



(31) 



T,. = f\4>) 

we can calculate the energy density for the electro-magnetic fields 

p_ = - < 0|T%|0 >= ^ / f 0<l' + ' (32) 



which can be estimated as 

4c-4 



a,; 



p,m -^HU-] , (33) 



where is the beginning of inflation. Once the energy density of the electro-magnetic fields 
Pem becomes comparable to the energy density of the inflaton namely, pem ^ Pcj, ^ MpHj, 
we need to consider the backreaction. More precisely, after the time a/aj ~ 10® we need to 
consider the backreaction. Usually, it is believed that once the electromagnetic fields becomes 
dominant, the backreaction spoils inflation. However, in the next section, we will see inflation 
continues opposed to the naive expectation. This is because the backreaction affects the 
dynamics of inflaton before the energy density of the electromagnetic fields dominates that 
of the inflaton. Thus, it is worth reexamining generation of magnetic fields from inflation 
with the backreaction. 



IV. MAGNETIC FIELDS FROM INFLATION WITH BACKREACTION 

We consider the potential V{(j)) = m^0^/2 in this section. Hence, we adopt the coupling 
function /(0) = e'^'^^''^*^? . Since we are going to look into the situation where the electro- 
magnetic field is not negligible, it is natural to consider an anisotropic spacetime with a 
coherent vector field. However, as the energy density of the vector field should be subdomi- 
nant during inflation, the anisotropy is negligible to the lowest order. Hence we ignore a in 
the basic equations and regard it perturbative quantity. 
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A. Inflation with backreaction 



The inflaton dynamics described by Eqs.(j4]) and ([7]) can be written as 

(34) 



1 

a — 



3M2 



-0 + -m + -e 



= -3d0 - + -^0e ^"^^ %i . (35) 

p 

When the effect of the vector field is comparable with that of the inflaton field as source 
terms in (1351) . we get the relation cp\t~'^'^'^ jM'^ ~ rr?. If we define the ratio of the 
energy density of the vector field pA = p\e~^'^'^^^p~'^"'/2 to that of the inflaton = m^0^/2 
as 

7? ^ ^ ^ , (36) 

we find the ratio becomes TZ ~ / ccf)^ when the above relation holds. Since inflation takes 
place typically at ~ (9(10)Mp, the ratio goes TZ ~ 10^^. Thus we find that the effect of 
the vector filed in flM|) is negligible even when it is comparable with that of the scalar field 
in dSS]). 

We notice that the above situation is not transient one but an attractor. Suppose that 
Pa is initially negligible, TZi ^ 10^^. In the conventional slow-roll inflationary phase ([9]), the 
relation e~'^'^l^^v oc t^"- holds as was shown in (fTOj) . Hence, the ratio TZ varies as 7?. oc ^^^'^-^)°' , 
As we now consider c > 1, increases rapidly during inflation and eventually reaches 
IZ ~ 10"^. Whereas, if the ratio is initially TZi ^ 10"^, the inflaton climbs up the potential 
due to the effect of the vector field in (!35|) . hence pA will decrease rapidly and again TZ ~ 10~^ 
will be realized soon. Thus irrespective of initial conditions, pA will track p^. 

From these arguments, the inflaton dynamics after tracking is governed by the modified 
slow-roll equations 

3d0 = -n?<\) + ^ • (38) 

p 

We refer to the phase governed by the above equations as the second inflationary phase, 
compared to the first one governed by the equations ([9]). Using above equations, we can 
deduce 

da ^ 
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This can be integrated as e-^<^'/*^'-^" = m'^M^{c - l)/c'^p\ [l + De"^^'^"^)"] ^ , where D is 
a constant of integration. This solution rapidly converges to 

Thus, we found pA becomes constant during the second inflationary phase [3] • In particular, 
from fl40l) . the relation 

/ = e^"^' oc (41) 
holds in this phase. Thus, when the backreaction becomes relevant, we effectively have 

Ceff = 1- 

We note that the backreaction makes the expansion anisotropic. In fact, as is done in a 
previous work [3], we can obtain a remarkable result from ([6]) as 

S cr 1 c — 1 , 

where e is a slow roll parameter and S/if is the anisotropic expansion rate normalized by 
the Hubble parameter. Typically, we have a small anisotropy S/iiT ~ 10"^. However, it is 
not impossible to detect this small number by PLANCK 



B. Generation of Magnetic fields 

In the situation where magnetic fields have a scale free spectrum c = 3/2, the formula 
( !2T|) tells us that the electric fields have a red spectrum. Hence, the largest scale has a 
dominant contribution to the energy density. We can assume the coherent electric fields 
with a definite direction dominate the energy density of the electromagnetic fields. Then, 
the result in the previous subsection is applicable. 

Now, we will take into account the backreaction. The point is as follows. Before the 
backreaction becomes important, we have the relation 

f^{^y\ (43, 

However, once the backreaction becomes important, we have an attractor behavior 

f^{±y\ (44) 
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That means the effective Ccs changes from Cefr = c to the critical value Ces = 1 due to 
backreaction. Since we have the formula (see Eq. fl33l) ) 

the transition point at occurs at pem ^ Pa ^ lO^^p^ ~ i.e. 



4c-4 



iq-'HjT'm; . (46) 



We are now in a position to calculate the power spectrum of magnetic fields. First, we 
consider the modes which exit the horizon before ai,. The superhorizon evolution of the 
mode function before aj, is given by 

MV) = T^(-T ' ' (47) 
fkV2k \akj 

where we should note fk is defined by = {ab/afY^~'^{ak/af)~^'^ from the continuity at at- 
Since the evolution after a?, becomes Wk oc a"^, we obtain the mode function after ah as 

..= ^(^V"(^V. (48) 



fkV2k \akj \ab 
From the formula (122|) . we obtain magnetic fields 



2 / -v \ 2c-3 / \ 2c-2 



Compared to the cases with no backreaction fl29l) . the amplitude is reduced by the factor 

)2c-2 / \ 2c-2 

- 10-'//-M„ (^) . (50) 

For the flat spectrum c = 3/2, we can deduce magnetic fields at the end of inflation as 

6B{X,.Vf) = lO-'Hy'M, (^) ^ . (51) 

\afj V27r 

Without backreaction, we anticipated lO^^^G for the scale invariant case c = 3/2. However, 
by taking into account the backreaction, we have a suppression factor a^/a/ in ( H9l) which 
is about 10~^^. Hence, we can expect at most 10~^^G on Mpc scales at present. For modes 
which exit the horizon after the transition time a^, by setting c = 1 in (H71) . we obtain 



3 
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where fk is now defined by fk = {ak/af) ^. Thus, we can calculate magnetic fields at the 
end of infiation as 

r2 / < X -1 



5b(Ap,?7/) 



Hf f X 



(53) 



V2n \Hy\ 

The resultant spectrum for c = 3/2 is schematically depicted in Fig.l. As is expected, we 
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-5* ^ I 



> k 



FIG. 1: The magnitude of magnetic fields is schematically depicted as a function of wavenumber k 
for the case c = 3/2. There exists a break at = abHj in the spectrum due to the backreaction. 
As can be seen, the amplitude of magnetic fields on Mpc scales gets a suppression by 10~24 due to 
this break. 

have a fiat spectrum on large scales before the backreaction becomes relevant. However, 
once the backreaction becomes important, the spectrum becomes blue. Thus, after taking 
into account the backreaction, we realized that primordial magnetic fields on large scales 
from infiation cannot be expected. 



V. CONCLUSION 



We studied cosmological magnetic fields in an infiationary scenario. We have explained 
how the backreaction works in magnetogenesis in the infiationary scenario. In usual analysis, 
it has been assumed that the backreaction spoils infiation. However, there is no justification 
for this assumption. Hence, we have incorporated the backreaction into the analysis and 
reanalyzed the production process of magnetic fields from infiation. As a consequence, we 
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have shown that inflation still continues after the backreaction becomes important. In spite 
of this fact, it turned out that creation of primordial magnetic flelds is still signiflcantly 
suppressed due to the backreaction. 

The main point is that the energy density of electric flelds have to grow in order to 
have sufficient amplitude for magnetic flelds, which causes the backreaction. What we have 
found is that the backreaction changes the dynamics of inflaton and, as a consequence, 
makes the genesis of magnetic flelds difficult. However, inflation is not spoiled by the back- 
reaction. Instead, an anisotropic inflationary universe has been created due to backreaction. 
Although magnetic flelds are negligible, since the expansion of the universe is anisotropic, 
other interesting phenomenologry such as the statistical anisotropy of primordial curvature 



perturbations can be expected 



9j. The possibility is now under investigation. 
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